Probing pairing gap in Fermi atoms by light scattering 
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We study stimulated scattering of polarized light in a two-component Fermi gas of atoms at 
zero temperature. Within the framework of Nambu-Gorkov formalism, we calculate the response 
function of superfluid gas taking into account the final state interactions. The dynamic structure 
factor deduced from the response function provides information about the pairing gap and the 
momentum distributions of atoms. Model calculations using local density approximation indicates 
that the pairing gap of trapped Fermi gas may be detectable by Bragg spectroscopy due to stimulated 
scattering. 
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I. INTRODUCTION 

Since the first realization of quantum degeneracy in an atomic Fermi gas by Jin's group |l| in 1999, cold Fermi 
atoms have been in focus of research interest in current physics. In a series of remarkable experiments, several 
groups 0, IE EL IE IE have demonstrated many new aspects of cold degenerate atomic Fermi gases. The ability 
to change interatomic interaction ranging from strong attraction to strong repulsion makes Fermi atoms the most 
favorable laboratory system for testing theoretical models in diverse fields. In particular, research with Fermi atoms 
has relevance in the field of superconductivity [E Il0| . The basic mechanism behind superconductivity and Fermi 
superfluidity is particle-particle pairing with an energy gap. Recently, two groups-Innsbruck and JILA 
have independently reported the measurement of pairing gap in Fermi atoms. Furthermore, two groups-Duke and 
Innsbruck [l3L IbH ] have measured collective oscillations which indicate the occurrence of fermionic superfluidity 01 
in atomic Fermi gas. The superfluid pairing is believed to occur near the crossover |T^ . IT7llT8| between the predicted 
BCS state of atoms and the Bose-Einstein condensation of molecules formed from Fermi atoms. Several groups have 
produced Bose-Einstein condensates (BEC) [l9j of molecules formed from degenerate Fermi atoms. Several other 
recent experimental (2(| and theoretical investigations [2lJ have revealed many intriguing aspects of cold Fermi gases. 

Several theoretical proposals [lEIl^ have been made for probing pairing gap. A method has been suggested to use 
resonant light | 22| to make an interface between normal and superfluid atoms. This has been recently implemented 
[Til |24| . albeit in the radio frequency domain. A number of authors |5E HE H3 have theoretically investigated 
Bogoliubov- Anderson (B A) mode [28|, I^E |3(| in fermionic atoms as a signature of superfluidity. B A mode constitutes 
a distinctive feature of superfluidity in neutral Fermi systems since it is associated with long wave Cooper-pair density 
fluctuations. 

Our purpose here is to study Bragg spectroscopy with off-resonant polarized lasers as a method for detecting the 
taring gap. Bragg spectroscopy has been used by Ketterle's group for measuring structure factor of an atomic BEC 
<H| . Bragg scattering in superfluid Fermi atoms has an analogy with Raman scattering in electronic superconductors 
|. In the next section we describe polarization-selective light scattering. We then discuss briefly the method of 
calculation of response function of superfluid Fermi gas using Green function techniques [^E 01 • We present the main 
results which suggest the possibility of detecting pairing gap by scattering of circularly polarized light. 



II. POLARIZATION-SELECTIVE LIGHT SCATTERING 



When two off-resonant laser beams with a small frequency difference are impinged on atoms, the scattering of one 
laser photon is stimulated by the other photon. In this process, one laser photon is annihilated and reappeared as a 
scattered photon propagating along the other laser beam. The magnitude of momentum transfer is q ~ 2k^ sin(6>/2), 
where 9 is the angle between the two beams and kr, is the momentum of a laser photon. To illustrate the main 
idea, we specifically consider trapped 6 Li Fermi atoms in their two lowest hyperfine spin states | g\) —\ 2Sx/2-,F = 
1/2, mp — 1/2) and | 52) =| 2Si/2,.F = 1/2, wif = —1/2). For simplicity, we consider that the number of atoms in 
each spin component is the same. An applied magnetic field tuned near the Feshbach resonance (~ 850 Gauss) results 
in splitting between the two spin states by ~ 75 MHz |35|. while the corresponding splitting between the excited 
states I e x ) =| 2P 3/2 ,F = 3/2, m F = -1/2) and | e 2 ) =j 2P 3/2 ,F = 3/2, m F = -3/2) is ~ 994 MHz 0. Let both 
the laser beams be cr_ polarized and tuned near the transition | 32) — H e 2 ). Then the transition between the states 
I gx) and | e-i) would be forbidden while the transition | gi) — >| ei) will be suppressed due to the large detuning 
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~ 900 MHz. This leads to a situation where the Bragg-scattered atoms remain in the same initial internal state | e 2 ). 
Similarly, atoms in state | 171) only would undergo Bragg scattering when two cr+ polarized lasers are tuned near the 
transition | g%) — »| 2P 3 / 2 ,-F = 3/2, rap = 3/2). Thus, we infer that in the presence of a high magnetic field, it is 
possible to scatter atoms selectively of either spin components only by using circularly polarized Bragg lasers. We 
assume that both the laser beams are cr_ polarized and tuned near the transition | gi) — >| €2). Under such conditions, 
considering a uniform gas of atoms, the effective laser-atom interaction Hamiltonian in electric-dipole approximation 
can be written as 

H T = m 7<r4(k + q)c CT (k)+H.c. (f) 

k,<j=l,2 

where c a (k) represents annihilation operator of an atom with momentum k in the internal state a. The subscript 
1(2) refer to the state | gi) (| 52)), ^ = (f2i + f^) /2, and ji = f2j/fi. Here fii denotes the two-photon Rabi frequency 
for the transitions | gi) — >| ei) — >| gi). For both the laser beams having <r_ polarization tuned near | 32) — H £2), we 
have Q2 » One can identify the operator p^\q) — ^ k c£.(k + q)cer(k) as the Fourier transform of the density 
operator. 



III. THE FORMALISM 



The scattering probability is given by the susceptibility 

X {^T-r') = -{T T [p^{r)p { 2 ] q {r')]). (2) 

where p^ — J2k a 'Y° a \+q o ak ,n T T is the complex time r ordering operator and (• • •) means thermal averaging. The 
dynamic structure factor is related to x by x( c li UJ n) as 

S(q,uj) = --[l + n B H]Im[ X (q,z = w+iO+)]. (3) 

TT 

This follows from generalized fluctuation-dissipation theorem. In order to treat collective excitations, it is essential to 
go beyond Hartree approximation and apply either a kinetic equation or a time-dependent Hartree-Fock equation or 
a random phase approximation 30] . The essential idea is to take into account the residual terms which arc neglected 
in the BCS approximation and thereby treat the off-diagonal matrix elements (vertex functions) of single-particle 
operators in a more accurate way |30l |36j . 

To study light scattering in Cooper-paired fermionic atoms, we apply Nambu-Gor'kov formalism [33ll3^| of super- 
conductivity 36]. Using the familiar Pauli matrices, the susceptibility can be expressed as 

X(q,w) = / ^0-Tr[ 7k G(k + )r(k + ,k_)G(k_)] (4) 
where the Green function has a matrix form as 

G{k) = k*-E? k +i6 ' (5) 
where E). = + and £/. = — p with e& = H 2 k 2 /(2m). Here To is a 2 x 2 unit matrix. The vertex equation is 

r(fc + ,M = 7 + i / ^T 3 G(fc;)r(fc;,fcL)G(fc^)r 3 u(k,k'), (6) 

where k± — k ± q/2 and k = (k, ko) is the energy-momentum 4- vector whose components are fc 3 = and k^ = iko- 
The bare vertex is a diagonal matrix: 7 = Diag.fyi, — 72]. Using Pauli matrices To and r 3 , this can be rewritten as 
1 = 7o(fc)To + 73(^)73, where 70 = [71 — 72] /2 and 7 3 = [71 + 72]/2. For unpolarized light in the absence of magnetic 
field, 71 =72- However, if the incident light is polarized, 71 7^ 72- In the specific case of <j_ polarization as discussed 
above, we have 7 3 ~ —70 — 72/2. If the potential V(k, k') is separable in k and k', then Eq. iJBJ is analytically 
solvable. We replace V(k, k') by the potential V — Airh 2 a s / '(2m), (m = m/2 being the reduced mass) expressed 
in terms of s-wave scattering length a s . The strong-coupling limit may be accessed by first renormalizing the BCS 
mean-filed interaction and then taking the limit a s — > ±00 as we will discuss later. The four-dimensional integrals 
of Eq. © can be performed following the established method of relativistic quantum electrodynamics as applied for 



3 



studying collective excitations in a superconductor |3j| . The detailed method of solution is discussed elsewhere [36 
We here present the final result 



X (q,c) = 2iV(0)7 2 <5)+2JV(0) 



where iV(0) is the density of states at the Fermi surface and 
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Here p q = ftq and v& is the velocity of the atoms with momentum k and 
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The symbol (X) implies averaging of a function X over the chemical potential surface: (X) = [N(0)] 1 J ci 3 k<5(efc)X. 
As to — > w + i0 + , /3 — > /9 + «0 + , we have the following analytic property of /(/3): 



sinhv//? 2 - 1 i7r/2 



, /9>1 



(10) 



IV. STRONG-COUPLING LIMIT 



To access the strong-coupling limit, the chemical potential [i and the gap should be obtained by solving the gap 
equation 



along with the equation 
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of the density of single component. Note that the eq. l|llf> is obtained by regularizing the zero-temperature BCS gap 
equation with a mean-field parameterized by two-body scattering length as done in Ref. [TtI ]. This approach fails to 
account for pairing fluctuation effects which are particularly significant near T c in strong-coupling regime. However, 
far below T c , the correction due to the pairing fluctuation is very small Based on this regularized mean-field 

approach and local density approximation (LDA), the zero-temperature density profiles |39j |. momentum distribution 
p0| and the finite temperature effects |4l| of superfluid trapped Fermi atoms have been recently studied . The two 
coupled eqs. Ijlll) and H12[) admit analytical solutions |4^. In the unitarity limit the solutions yield A ~ 1.16m, 
(j, = (1 + j3)ep 1 where (3 — —0.41 [ioll42]| is a constant. Many other recent theoretical [43T,l44| and experimental 0,^3 
studies have established the universality of Fermi gas in the unitarity limit. 



V. LEADING APPROXIMATIONS 

With the use of Eqs. (|7I10|) in Eq. J3J), the dynamic structure factor in the leading approximation in terms of /3 _1 
can be written as 



S(q,u)=N(0) 
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(3 > 1 



(13) 



This is also obtainable from the BCS- Bogoliubov mean-filed treatment as shown in [38|. Although this leading 
approximation takes into account excitations in the particle-hole continuum, it fails to account for the in-gap collective 
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FIG. 1: Dimensionless dynamic structure factor S(S, q)/N(0) of a uniform superfluid Fermi gas is plotted as a function of 
dimensionless energy transfer uj/eo (eo is the Fermi energy) for different values of the scattering length a a = 2.76k p 1 (solid), 
a a — S&dkp 1 (dotted), a s — 5A7kp 1 (dashed) for a fixed momentum transfer q — 0.8/cf- The dash-dotted curve is plotted for 
a s = 2.76kp 1 and q = OAkp. The inset shows the variation of the gap A and the chemical potential n as a function of a a . 
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FIG. 2: Same as in Fig. (|5J but for a trapped superfluid Fermi gas for a fixed momentum transfer q — 0.8&F. Here eo is the 
Fermi energy at the trap center 



modes. It can be verified [2S| that in the regime of large momentum and energy transfer ((3 » 1), the dynamic 
structure factor of Eq. approximately satisfies the f-sum rule 

o;5(q,a;)da;~^, £q»l (14) 
2m 

where N is the total number of particles. 

Now let us consider the case < f3 « 1, that is v^.p^ < to « 2 A. In this case, the second term in the coefficient 
of 73 in Eq. J7J dominates over all other terms. This term leads to BA mode appearing as a pole of \. BA mode 
restores the continuous symmetry that is broken by BCS ground state. In the limit q — > and u) — * 0, / ~ 1 and 
hence the pole is 

1 . . 

WBA = -j=v F p q . (15) 

In the low momentum and low energy limit (0 < /3 « 1) the dynamic structure factor can be obtained by linearizing 
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the denominator of the second term in Eq. J7J around the BA mode. By approximating / ~ 1, we then obtain 

S(q, u) = N(0)^^—5(u - lu ba ). (16) 
Zuj ba 

With 73 — > 1, this satisfies the /—sum rule. BA mode is well defined in the low momentum regime, i.e., for £q = 
vpPq/ (2A) « 1. For large momentum, it becomes ill defined due to Landau damping. Ohashi and Griffin [26| have 
provided a detailed theoretical treatment of this mode in the BCS-BEC crossover in Fermi atoms. Minguzi et. al. 
|27| have found that this mode appears as a prominent asymmetric peak in the spectrum of density fluctuation at a 
very low momentum and energy. 



VI. RESULTS AND DISCUSSIONS 



Figure 1 and 2 show S(lu, q) as a function of u> for a uniform and trapped gas, respectively, for different values of 
a s . In the case of trapped gas, we use LDA with local chemical potential //(r) determined from equation of state of 
interacting Fermi atoms in a harmonic trap. When a s is large, the behavior of S(5, q) is quite different from that of 
normal as well as weak-coupling BCS superfluid. This can be attributed to the occurrence of large gap for large a s . 
In contrast to the case of a uniform superfluid [32|, S(S, q) for a superfluid trapped Fermi gas has a structure below 
2A(0), where A(0) is the gap at the trap center. As the energy transfer decreases below 2A(0), the slope of S(6,q) 
gradually reduces. Particularly distinguishing feature of S(S, q) of a superfluid compared to normal fluid is gradual 
shift of the peak as a s or A increases. The quasiparticle excitations occur only when 2A(x) < uj. This implies that, 
when to is less than 2A(0), the atoms at the central region of the trap can not contribute to quasiparticle response. 



VII. CONCLUSIONS 

Order of magnitude analysis of Ref. |46j suggests that, with large momentum transfer, it may be possible to 
distinguish the scattered atoms in time of flight images. A comparison of images with and without Bragg pulses 
may reveal information about the momentum and density distribution of the scattered atoms. Furthermore, the 
polarization-selective Bragg spectroscopy may lead to better precision in time-of-flight spin-selective measurements 
since they will be in the same spin component. It is possible to select counter propagating scattered atoms using 
three or four beam scattering configurations 46]. One can then explore the possibility of measuring the correlation 
of two scattered atoms with opposite momentum by the technique as used in recent studies |47l l48j . 
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